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The ground state properties of the high spin Heisenberg chains with alternating single site 
anisotropy are investigated by means of the numerical exact daigonaization and DMRG method. 
It is found that the ferrimagnetic state appears between the Haldane phase and period doubled 
Neel phase for the integer spin chains. On the other hand, the transition from the Tomonaga- 
Luttinger liquid state into the ferrimagnetic state takes place for the half-odd-integer spin 
chains. In the ferrimagnetic phase, the spontaneous magnetization varies continuously with the 
modulation amplitude of the single site anisotropy. Eventually, the magnetization is locked to 
fractional values of the saturated magnetization. These fractional values satisfy the Oshikawa- 
Yamanaka- Affleck condition. The local spin profile is calculated to reveal the physical nature of 
each state. In contrast to the case of frustration induced ferrimagnetism, no incommensurate 
magnetic superstructure is found. 

KEYWORDS: high spin Heisenberg chain, single-site anisotropy, spatial modulation, Haldane phase, ferri- 
magnetic phase, period-doubled Neel phase, hidden antiferromagnetic order 
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1. Introduction 

Among various exotic ground states in quantum mag- 
netism, the Haldane state in the integer spin antiferro- 
magnetic Heisenberg chain 1 has been most extensively 
studied both experimentally and theoretically. This state 
is characterized by the hidden antiferromagnetic string 
order accompanied by the Zi x Z2 symmetry breakdown 
in spite of the presence of the energy gap and exponen- 
tial decay of the spin-spin correlation function. The easy 
plane single-site anisotropy D{> 0) destroys the Haldane 
ground state leading to the large- D state with finite en- 
ergy gap and exponential decay of the spin-spin corre- 
lation function without specific order. On the contrary, 
the easy axis single-site anisotropy (D < 0) drives the 
Haldane state into the Neel state. 2-4 

On the other hand, the ground state of the odd spin 
Heisenberg chain is the Tomonaga-Luttinger liquid state. 
Due to its critical nature, the ground state is driven to 
the Neel ordered state for infinitesimal negative D while 
the Luttinger liquid state is stable against positive D. 5 

In this context, it is an interesting issue to investi- 
gate how the ground states of the quantum spin chains 
are modified if the easy-axis and easy-plane single-site 
anisotropy coexist in a single chain. In the previous 
work, 6 the present author and Chen investigated the 
5=1 chain with alternating single-site anisotropy and 
found that the period doubled Neel phase with j J. 0) 
structure is realized for strong alternation amplitude, 
although the Haldane phase is stable for weak alter- 
nation. The physical origin of this type of Neel order 
is interpreted as a 'pinning' of the string order. In the 
present work, we further explore this problem for the 
cases 5 > 1. We find not only the || J, 0) ground 
state but also the ferrimagnetic ground states with quan- 
tized and unquantized spontaneous magnetization for 
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intermediate strength of I?- alternation. These quanti- 
zated values of magnetization also satisfy the Oshikawa- 
Yamanaka- Affleck condition 7 well-known for the magne- 
tization plateau in the magnetic field. 

This paper is organized as follows. In the next section, 
the model Hamiltonian is presented and the two possible 
senarios which leads to different ground states are ex- 
plained. In §3, the numerical results for the spontaneous 
magnetization and the local spin profile are presented to 
reveal the physical nature of each state. The last section 
is devoted to summary and discussion. 

2. Model Hamiltonian 

We investigate the ground state of the Heisenberg 
chains with alternating single site anisotropy whose 
Hamiltonian is given by, 



n = JSlS > 
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(J > 0,SD > 0). 



(1) 



where Si is the spin- 5 operator on the z-th site. 

In ref. 6, it is found that the period doubled Neel 
phase with || ! 0) structure is realized for large enough 
8D for 5 = 1, although the Haldane phase is stable for 
small SD. The mechanism to stabilize this period dou- 
bled Neel phase can be understood along the following 
senario (senario I). In the absence of the Z?-terms, the 
5 = 1 ground state has a hidden string order which 
implies that the spins with |±1) are arranged antifer- 
romagnetically if the sites with |0) are skipped. 2,3 The 
position of the sites with |±1) and |0) strongly fluctuate 
quantum mechanically and this antiferromagnetic order 
remains hidden because it is impossible to observe the 
correlation between only the sites with |±1) experimen- 
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Fig. 1. The classical spin configuration in the ferrimagnetic state 




Fig. 2. The spontaneous magnetization in the classical limit. Here 
and in the following figures energy scale is set as J = 1. 



tally. In the presence of strong <5D-terms, only the states 
consistent with the constraint set by these SD-teims sur- 
vive among all states with hidden order. For SD » J, 
the odd-th site must be in the state |0) and the even-th 
sites |±1). To be compatible with the string order, spins 
must be arranged as 1 1 J. 0) . Thus the strong SD-tevm 
can select the spin states among those with hidden order 
to realize the explicit period doubled Neel order. 

On the other hand, another senario (senario II) is 
possible from classical point of view, although this is 
not realized in the S = 1 case. Let us consider the 
classical limit where each spin can be regarded as a 
classical unit vector. In the ground state, the spins on 
the easy axis sites are fully polarized along ^-direction 
but those on the easy-plane sites are tilted by an an- 
gle 9 = cos -1 (J /SD) from the — z-dircction as shown in 
Fig. 1. Therefore this senario leads to the noncollinear 
ferrimagnetic ground state with the spontanenous mag- 
netization M = M s (l - (J/5D))/2 as plotted in Fig. 2. 

In what follows, we show that either of these two senar- 
ios can be realized in S > 1 chains depending on the 
values of SD and S based on the numerical exact di- 
agonaliztion (NED) and density matrix renormalization 
group (DMRG) calculation. 

3. Numerical Results 

3.1 Spontaneous Magnetization 

To identify the ferrimagnetic regime expected in 
senario II, the ground state spontaneous magnetization 
is calculated by NED with periodic boundary condition 
and by DMRG with open boundary condition for various 
values of S and SD. The maximum chain length for the 
NED is TV = 12 for S = 3/2, S = 2 and S = 5/2, while it 
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Fig. 3. The spontaneous magnetization for (a) 5 = 2 with N = 
12 (NED) and N = 64(DMRG), and (b) S = 3 with N = 8 
(NED) and N = 32 (DMRG) plotted against &D. The dotted 
lines are the classical spontaneous magnetization. 



is N = 8 for 5" = 3. For the DMRG calculation with open 
boundary condition, appropriate end spins are added to 
reduce the boundary effects. 

The results for the integer spin cases are presented in 
Fig. 3 for 5* = 2 and 3 and those for the half-odd-integer 
spin cases are presented in Fig. 4 for S — 3/2 and 5/2. 
In contrast to the case of S = 1, it is found that the 
ferrimagnetic phase always appear for S > 3/2 above 
the critical value SD c i. For < SD < SD c i, the energy 
gap decreases monotonously with SD until it vanishes at 
SD = SD cl in all cases studied (3/2 < S < 3). There- 
fore, we may safely conclude that the ground state is the 
Haldane phase or the Tomonaga-Luttinger liquid phase 
according as 5* = integer or half-odd-unteger. 

For integer S, the spontaneous magnetization vanishes 
for large enough SD. From the local spin profile (S 1 ?) 
which will be presented in the next section, this state 
turns out to be the period-doubled |f J. 0)-type Neel 
state expected in senario I. On the other hand, the ferri- 
magnetic state remains stable for arbitrarily large SD for 
half-odd-intcger S because the ground state of the easy 
plane site is the doublet S z — ±1/2 which can sustains 
the ferrimagnetic order even for large SD. 

It should be also noted that the spontaneous magne- 
tization in the ferrimagnetic phase is not restricted to 
simple fractions of the saturation magnetization M s (= 
NS) as in the usual quantum ferrimagnets 8,9 but varies 
continuously with SD in accordance with the classical 
intuition predicting the noncollinear ferrimagnetism in 
senario II. 

Within a appropriate range of SD, however, the spon- 
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Fig. 4. The spontaneous magnetization for (a) S = 3/2 with TV = 
12 (NED) and N = 64(DMRG) and (b) S = 5/2 with N = 12 
(NED) and N = 32(DMRG) plotted against SD. The dotted 
lines are the classical spontaneous magnetization. 



taneous magnetization is locked to a simple fraction of 
M s reflecting the quantum nature of the present model. 
In the DMRG calculation, these quantized value of mag- 
netization slightly deviates from the simple fraction of 
M s due to the boundary spins. Correspondingly, the 
spontaneous magnetization is slightly rescaled so that 
the main quantized value exactly equals the simple frac- 
tion of M s in Figs. 3 and 4. 

In all cases, these quantized values of the magnetiza- 
tion satisfy the condition 



p(S -m) = q 



(2) 



where p is the size of the unit cell, q is an integer and 
m is the magnetization per site (m = M/N = MS/M S ). 
In the present model, p is equal to 2. This condition is 
identical to that proposed by Oshikawa, Yamanaka and 
Affleck 7 for the magnetization plateau in the magnetic 
field. However, their proof is not restricted to the mag- 
netic field induced magnetization but also applies for the 
spontaneous magnetization in the ferrimagnetic phase. 
If the condition (2) is satisfied, it is allowed to have a 
finite energy gap to the excited state with different mag- 
netization. This implies the stability of the ground state 
against the variation of SD which leads to the 'plateau' 
behavior. The <5D-depcndcnce of the energy gap on the 
plateau state calculated by the DMRG method for 5 = 2 
chains is shown in Fig. 5. It is clear that the energy gap 
is finite on the plateau region 1.91 < SD < 3.26. 

In addition, the maximum possible values of the 
ground state spontaneous magnetization is bounded from 
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Fig. 5. The energy gap of S = 2 chain on the plateau state with 
magnetization M p = M s /4 for N = 72. The filled (open) symbols 
are the gap to the state with magnetization M = M p + l(M p — 1) 



above due to the nature of the present model. To max- 
imize the spontaneous magnetization, the spins on the 
easy axis sites must have Sf = 5 and those on the easy 
plane site must have negative value due to the antifcr- 
romagnetic interaction with the neighbouring polarized 
spins. It should be noted that the senario I leading to the 
period doubled Neel order becomes effective if Sf = on 
the easy plane site. Therefore the spontaneous magneti- 
zation per site m must satisfy < m — MS/M S < 5/2 
in the ferrimagnetic phase. This implies 25 > q > 5. 
In the half-odd-integer case, the smallest possible value 
of q is 5 + 1/2. This gives the magnetization per site 
m = (5 - l/2)/2 which yields M/M s = (5 - l/2)/25. 
On the other hand, in the integer spin case, the smallest 
possible value of q is 5 + 1. In this case, m is equal to 
(5 - l)/2 which yields M/M s = (5 - l)/25. It should be 
noted this value vanishes for 5 = 1. This explains why 
the quantized ferrimagnetic phase does not appear for 
5=1 case. 

Actually, prominent 'plateaus' are observed only for 
the smallest possible value of q. This is due to the fact 
that the condition for the gap generation on the com- 
pactification radius of the underlying Gaussian model 
becomes increasingly severer with the increase of q. 7 As 
a secondary plateau with larger q, we only find a small 
plateau at M = M s /5 for 5 = 5/2 which corresponds to 
q = 4(= 5 + 3/2) as shown in Fig. 6 within the 5 values 
studied so far, although this plateau is almost unvisible 
for small sized systems shown in Fig. 4(b). 

3.2 Local Magnetization Profile 

The local magnetization profile (5?) calculated by the 
DMRG method is presented for each phase of 5 = 2 
chains in Fig. 7. Below the plateau, the easy plane spins 
are almost in the state Sf = —1 while the magnetiza- 
tion of the easy axis spins increases from 1 to 2 as SD 
approaches the lower end of the plateau region. On the 
plateau, the easy axis spins are almost in the state Sf = 2 
and the easy plane spins are in the state Sf = —I lead- 
ing to the quantized value of spontaneous magnetization 
for the smallest possible value of q described in the pre- 
ceding section. Above the plateau, the easy axis spins 
are almost in the state Sf = 2 and the increase in total 
spontaneous magnetization is due to the decrease in the 
polarization of the easy plane spins. 
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Fig. 6. The small 'plateau' at spontaneous magnetization M = 
M s 1 5 for 5 = 5/2 plotted against 5D . The system size is W = 96 
(dotted line) and 192 (solid line). Rescaling factor is slightly 
different from that of Fig. 4 to fix this small plateau precisely to 
M = M s /5. 

The behavior of the local magnetization profile in the 
noncollinear ferrimagnetic phase is in contrast to the sim- 
ilar noncollinear ferrimagnetic phase in the frustrated 
spin chains investigated in refs. 10 and 11 in which the 
incommensurate superstructure is observed. This sug- 
gests that the incommensurate superstructures observed 
in these literatures are essentially due to frustration. 

In these ferrimagnetic phases, the correlation between 
the easy axis spins are ferromagnetic. At SD = 5D C 2, 
however, the easy plane spins turns into the state with 
Sf = and the correlation between the easy axis spins 
turns into antiferromagnetic. In this case, the magnetion 
profile clearly shows the |t J. 0) structure as shown in 
Fig. 7(d). For the calculation of the local spin profile in 
this phase, we have applied a tiny symmetry breaking 
field with period 4, because otherwise the true ground 
state of the finite size system is the linear combination 
of 1 1 | 0} and | J, \ 0)-type states and no net local mag- 
netization is expected. Actually, the DMRG calculation 
is often trapped to the states with domain walls in the 
absence of the symmetry breaking field. The value of the 
symmetry breaking field ranged from 0.005 to 0.02 and 
the results turned out to be almost indistinguishable on 
the scale of the Fig. 7(d). The physical origin of this 
magnetic structure is understood in the same way as the 
S = 1 case 6 following the first senario described in §2. 

4. Summary and Discussion 

The ground state properties of the high spin Hcisen- 
berg chains with alternating single site anisotropy are 
investigated by means of the numerical exact daigonaiza- 
tion and DMRG method. It is found that the ferri- 
magnetic state appears between the Haldane phase and 
period doubled Neel phase for the integer spin chains. 
On the other hand, the transition from the Tomonaga- 
Luttingcr liquid state into the ferrimagnetic state takes 
place for the half-odd- integer spin chains. In the ferri- 
magnetic phase, the spontaneous magnetization varies 
continuously with the modulation amplitude of the sin- 
gle site anisotropy in accordance with the classical in- 
tuition. Eventually, however, the magnetization is locked 
to fractional values of the saturated magnetization which 




satisfies the Oshikawa-Yamanaka- Affleck condition. 

The local spin profile is calculated to reveal the phys- 
ical nature of each state. In contrast to the case of frus- 
tration induced fcrrimagnetism, 10, 11 no incommensurate 
superstructure is found. We thus expect that the incom- 
mensurate superstructures found is these literatures are 
essentially due to the interplay of quantum effect and 
frustration. 

The similar mechanism should also work in 2 and 3 di- 
mensions, although the Haldane or Tomonaga Luttingcr 
liquid phase would be replaced by the long range ordered 
Neel- type state. However, even in the large SD limit, 
the ground state is not trivial due to the frustration in 
the effective interaction among the easy-axis spins. The 
investigation of these higher dimensional models are in 
progress. 
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For the experimental realization of the present mech- 
anism it is necessary to synthesize the compound of easy 
axis magnetic ions and easy plane ones. Considering a 
variety of phases expected for the present model, this is 
a challenging attempt for the experimentalists. Recently, 
various single chain molecular magnets with consider- 
able strength of single site anisotropy have been synthe- 
sized. 12,13 Although the materials with alternating sign 
D-terms with uniform S are not yet reported, these se- 
ries of materials can be a good candidate to observe the 
phenomena proposed in the present work. 

The computation in this work has been done using 
the facilities of the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo and the Infor- 
mation Processing Center of Saitama University. The di- 
agonalization program is based on the TITPACK ver.2 
coded by H. Nishimori. This works is supported by the 
Grant-in-Aid for Scientific Research from the Ministry 
of Education, Culture, Sports, Science and Technology, 
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